Scaling of interfaces in brittle fracture and perfect plasticity 



O 

o 
o 

(N 



00 

o 

B 



i 

-a 
c 

o 
o 



> 

m 

m 
o 
o 
o 

-I— > 



i 

c 

o 
o 



X 



Eira T. Seppala, Vilho I. Raisanen, and Mikko J. Alava 
Helsinki University of Technology, Laboratory of Physics, 
P.O.Box 1100, FIN-02015 HUT, Finland 
(February 1, 2008) 

The roughness properties of two-dimensional fracture surfaces as created by the slow failure of 
random fuse networks are considered and compared to yield surfaces of perfect plasticity with similar 
disorder. By studying systems up to a linear size L — 350 it is found that in the cases studied 
the fracture surfaces exhibit self-affine scaling with a roughness exponent close to 2/3, which is 
asymptotically exactly true for plasticity though finite-size effects are evident for both. The overlap 
of yield or minimum energy and fracture surfaces with exactly the same disorder configuration is 
shown to be a decreasing function of the system size and to be of a rather large magnitude for all 
cases studied. The typical "overlap cluster" length between pairs of such interfaces converges to a 
constant with L increasing. 
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I. INTRODUCTION 



Roughness of fracture surfaces (FS) is a currently top- 
ical problem that has opened up surprising connections 
between engineering and weakly understood questions 
of statistical physics. The simple question why and 
how a crack surface becomes rough dodges easy answers 
since there is a multitude of experimental facts and ways 
for cracks to develop or propagate. One simplification, 
adopted in this work is to neglect cases in which the 
prevalent feature is the propagation of "fast" cracks in 
the favor of slow, adiabatic crack formation. The ques- 
tions we address here are related to how disorder affects 
crack surfaces and how interfaces born out of different 
load-elongation responses are related. Disorder is present 
in materials at all length scales in the form of atomic im- 
purities, dislocations, grain boundaries and so forth. 

No generally accepted picture exists yet of how slow 
cracks are formed and how it relates to crack interfaces 
tO. In three dimensions there are indications that the 
cracks become above a certain intermediate length scale 
self-affine so that the roughness exponent £ is close to 
0.8. Moreover the physics of crack advancement indi- 
cates that the generic features of phase transitions of 
driven lines (crack fronts in three dimensions) become 
relevant 0]. Quantitative agreement is however missing. 
For slow fracture in two dimensions (2D) and at small 
length scales in 3D the interface scaling may be differ- 
ent in that the exponents are close to those of the min- 
imum energy (ME) interface. These are the same as for 
the Random Exchange Ising Model (REIM) domain walls 
at zero temperature, and have therefore the exact value 
C = 2/3 in 2D g and the approximate value 0.41 ± 0.01 
in 3D 0-^]. The physics involved is simple: the crack 
minimizes up to the pertinent length scale the surface 
energy E given by 



E = J d*- 1 * [r{Vz(x)} 2 + l/ r {x,z(x)}] , 



(1) 



where the integral counts in two contributions. One 
arises from a surface stiffness (proportional to T) due to 
the deviations V-z(x) from a straight crack, and a second 
from a random disorder potential with a two-point cor- 
relator (V r (x' , z')V r (x, z)) where the disorder average is 
implied and (x', z'),(x, z) denote two locations inside the 
medium. The exponents quoted above are true in case 
the disorder has point-like correlations. The fluctuations 
of the potential would in an experiment correspond to 
a varying failure threshold or elastic modulus etc., de- 
pending on the circumstances. The closeness of the nu- 
merical values of the roughness exponents gives rise to 
the intriguing question as to why should a slow fracture 
resemble a global optimization like ground state domain 
walls? The connection is suggested by the fact that in 
arbitrary dimensions, lattice models that describe scalar 
perfect plasticity can be exactly mapped to the REIM 
domain wall problems. For brittle fracture or vectorial 
failure problems in general the correspondence is not ob- 
vious. Two-dimensional failure is special in that there is 
some experimental evidence of the crack roughness scal- 
ing with the domain wall in REIM, i.e. the so called 
directed polymer (DP), global roughness exponent [Q. 
This connection between global optimization and frac- 
ture surfaces has also been made in 2D simulations of 
brittle failure J8|||. In 3D, it is still unclear whether even 
just numerical models show such a universality [p|-|Tl"|]. 

In this paper we investigate in two dimensions the 
scaling properties of slow fracture surfaces and compare 
them to minimum energy surfaces with similar a priori 
disorder. We perform numerical simulations of the ran- 
dom fuse network (RFN) model, which has been stud- 
ied extensively as a model of brittle failure of disordered 
materials [p^|-^4[. As a warm up, we consider exten- 
sive system properties such as fracture stress, fracture 
strain, and damage. For the main case studied here, 
the dilution-type disorder, these are found to be in good 
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agreement with the critical defect -type arguments pro- 
posed by Duxbury and co-workers |13|] that imply loga- 
rithmic scaling with system size. Note that when com- 
paring with "reality" this kind of models contain two as- 
sumptions: first, the stress relaxation is supposed to be 
much faster than the stress rate (an adiabatic failure). 
Second, one assumes that the energy released by local 
crack formation is dissipated with no effect on the crack 
propagation. 

The paper starts with a short description of the nu- 
merical methods used and the dynamics of the adiabatic 
crack formation in Section [n| Section III discusses the 
strength properties of random fuse networks as a function 
of system size L. There are a number of ways to char- 
acterize a posteriori a self-affine interface. This is the 
main theme of Section [V, the topic of showing that 2D 
brittle fracture interfaces have the DP-type scaling. We 
demonstrate how both the so-called local width and the 
statistical properties of ensembles of interfaces indicate a 
similar kind of self-affine scaling. The scaling exponent is 
seen to be close to the DP one, £ = 2/3. The section also 
contains numerical data for varying disorder strength, 
and in particular compares perfect plasticity and brittle 
fracture by measuring the overlap of the associated in- 
terfaces starting with the same disorder configurations. 
This would be particularly relevant should it be so that 
the fracture and yield interfaces use the same "valley" in 
the landscape of the energies or thresholds. The paper is 
finished with a discussion in Section 0. 



II. CREATING THE INTERFACES 

A. Numerical models 

Random fuse networks are electrical analogues of elas- 
ticity and failure with disorder incorporated. One usually 
sets to mimic a tensile test implying that the extensive 
thermodynamic parameters become V ex t and I ex t, exter- 
nal voltage and current, respectively. These correspond 
to displacement and force in a real experiment. To study 
brittle failure one defines the elements that connect two 
nodes on an original lattice as fuses. These have a linear 
voltage-current relationship until a breakdown current if,, 
see Fig. A second choice would correspond to perfect 
plasticity if one made the fuses such that the local current 
becomes irreversibly constant at i y and stays so unless 
the local voltage is reduced, in which case the conduc- 
tivity becomes the original one and there is a permanent 
yield strain. 

In the following we use two different numerical tech- 
niques to study both brittle and perfectly plastic RFN's. 
Brittle failure is studied with standard adiabatic fracture 
iterations. These consist of solving the current balance in 
the system from Kirchoff 's and Ohm's laws and breaking 
after each iteration the most strained fuse (the criterion 



is min(ij/ J Cj j), where ij is the local current in each of 
the fuses and J CJ is the local threshold). The currents 
and voltages are found by solving the linear system of 
currents by the conjugate-gradient method. 

For perfect plasticity we use a mapping to minimum 
energy interfaces, i.e. random exchange Ising domain 
walls in their ground state, where exchange constants 
Jij between nearest neighboring spins are random but 
non-negative. In some cases we have exactly the same 
quenched disorder (equal thresholds for failure it and 
yielding i y = 2Jij for each fuse) as for brittle failure, and 
in the following the threshold for a fuse in both cases is 
denoted by J c . The simulations are done using combi- 
natorial optimization: finding the yield path (in 2D) is 
equivalent to the minimum cut-maximum flow problem 
of network flows |l5| ] that minimizes X^mter/ace *V This 
technique is more convenient than the transfer matrix 
methods in that there are no restrictions for the shape of 
the optimal path as overhangs and arbitrary transverse 
steps arc included in a natural fashion. 

The typical choice for introducing disorder to a RFN 
is to pick the failure currents J c from a prescribed prob- 
ability distribution P(J C ). The important issue is the 
behavior of P for J c ~ and for J c — > oo; the tails 
of the distribution are known to have strong effects on 
the strength properties and damage accumulation in the 
case of brittle fracture. For perfect plasticity or directed 
polymers the case is much simpler in that for one di- 
mensional interfaces in (1 + 1) dimensional systems such 
point-like disorder is in the renormalization group sense 
always relevant. Thus one expects always the same scal- 
ing properties in terms of interface roughness and samplc- 
to-sample interface energy fluctuations, these correspond 
to yield stress fluctuations in plasticity. The amplitudes 
are however non-universal and thus will depend on the 
exact form of P. 

In the following we study as typical examples the cases 
where P(J C ) is a flat distribution (P(J C ) = 1/(26 J) for 
Jo — 5 J < J c < Jq + SJ) and where P corresponds to 
"dilution disorder". That is P(J C ) = pS(J c - 1) + (1 - 
p) 8(J C ). The fraction of fuses that remain for infinitesi- 
mal currents with dilution is denoted with p, which has 
a value p — 0.8 unless otherwise mentioned, like for the 
uniform distribution case 5J/Jq = 1. The systems are 
chosen so that the direction of macroscopic current flow 
is aligned in the (10) orientation of the square lattice hav- 
ing periodic boundaries in the perpendicular direction. 
The systems are isomorphic, i.e. L x = L z , and their 
sizes range from L? = 10 2 to 350 2 for brittle failure and 
to 1000 2 for perfect plasticity, respectively. The mean 
positions of the surfaces are not fixed, hence they may 
sit anywhere in the system. The interfaces are defined 
in the usual way so that in the case of overhangs, the so 
called solid-on-solid approximation is used, i.e. the inter- 
face is found by taking the furthest value of the interface 
with respect to a fixed end of the network. The num- 
ber of realizations N over which the disorder-averaging 
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is performed is limited by the CPU-time for the simu- 
lations of the brittle fracture. In the case of plasticity, 
the technique used leads to a roughly linear scaling with 
respect the number of fuses in a system regardless of the 
threshold distribution. The number of different random 
realizations is shown in Table | for the cases in which 
exactly the same random networks are studied for both 
the brittle failure and the perfect plasticity. If only the 
ME-surfaces are studied N = 200 - 5000. 



how the field of local stresses will develop and normalize 
the local thresholds %j in the failure criterion. Consid- 
ered in this light, it is sensible that the brittle fracture 
surfaces are "blocking paths" , too. Yet the question re- 
mains whether the interfaces are still in the same uni- 
versality class: if the correlations in the renormalized 
disorder become different enough from point-like corre- 
lations the interface scaling properties will change. E.g. 
columnar correlations (U r (x, z) constant along Xi or z) 
would be relevant in this respect. 



B. Formation of the interfaces 



III. SCALING OF FRACTURE 



One should note that if fracture surfaces have non- 
trivial geometric scaling properties and in particular re- 
semble directed polymers (in 2D) this opens up several 
further questions: whether the outcome is independent of 
disorder strength, whether the disorder is always relevant 
for 2D minimum energy surfaces, and how the surface 
roughness relates to other typical quantities. The stan- 
dard way of iterating fracture in fuse networks, whether 
perfectly plastic or brittle, is based on extremal dynam- 
ics. The condition min(ij/ J c ,j) for the failure of the next 
element contains two effects: the disorder through the 
threshold and the local current that depends on the en- 
vironment of the fuse. 

For perfect plasticity, the information necessary for 
finding the final yield surface is contained in the initial 
field J c ,j due to the monotonicity property noticed by 



Roux and Hansen |16|. Even if one simulates the devel- 
opment of the system as a serie of fuse network problems 
(the tangent problem), the local current never decreases 
in a yielding process. Thus the final yield surface equals 
a blocking configuration that can be calculated from the 
original thresholds. This is related to the fact that the 
surface is much faster to compute than the whole pro- 
cess by considering it as an optimization problem for the 
interface: the history of the whole process or system in- 
volves much more information. 

For brittle fracture the monotonicity property is not 
true and thus no direct mapping exists between the initial 
disorder and a quantity to be minimized. The mapping of 
perfect plasticity to fuse networks makes it on the other 
hand clear that the difference between the processes is 
smaller than it would seem at the first glance. This is 
because in the tangent algorithm one has to solve a serie 
of adiabatic failure problems with the local yield thresh- 
olds i v renormalized by subtracting the current already 
passing through the fuse. Nonetheless each failure itera- 
tion is affected by stress-enhancement effects exactly as 
in a failure problem with the same fuses still intact. For 
brittle failure, the implication of the stress-enhancements 
during the failure process is that in order to obtain a min- 
imum energy surface (as defined by Eq. (|l|)) the original 
disorder ibj has to be renormalized. That is, the thresh- 
olds or missing fuses contain frozen-in information about 



A fracture can be contrasted with perfect plasticity 
also by looking at extensive thermodynamic quantities. 
The standard ones to consider are the damage rif,, the 
number of fuses broken in total, and the failure current lb 
and voltage Vb as computed from the maximum current 
of the TV-curve. For "truly" brittle failure this definition 
of Vb agrees with with that defined as the end point of the 
TV-curve. In the failure of brittle fuse networks there is 
considerable evidence for the relevance of critical defect 
-type effects. That is, the defect with the largest current 
enhancement will dictate the scaling of the current and 
voltage. For yield surfaces one would have I y = E ~ L, 
AE r~j L e where 8 = 2( — 1 and ( is the roughness expo- 
nent. Thus the critical strength quantities, without the 
renormalization discussed above, are supposed to have 
different scaling behavior in plasticity and fracture. 

Fig. H shows the scaling of damage and the strength 
quantities for dilution-type disorder, p = 0.8. The lines 
in the figure have been found with least-square-fits to 
data using the assumptions of linear scaling for n& and 
for the other two quantities the scaling Vb, lb ~ L/\/\n~L, 
which comes from the extreme value statistics, i.e. Gum- 
bel distribution, studied in the fracture case by Duxbury 
and coauthors |U| . It is seen that the scaling of the num- 
ber of broken fuses is asymptotically very close to a linear 
one. This means that the system still breaks in a brittle 
mode for p = 0.8. For Vj, and lb the scaling in the whole 
regime follows beautifully the Lj\J\nL scaling. Notice 
that the surface energy of yield surfaces is in principle a 
lower limit for ra& and that both scale linearly. One sees 
that the energy of yield surfaces or the lower limit for rib 
is lower by a constant factor than rib of fracture surfaces. 
Similar behavior is visible in the roughness values of the 
fracture and yield surfaces studied in the next section. 



IV. SCALING OF INTERFACES 



A. Global and local interface width 



There are several ways to characterize the scaling prop- 
erties of interfaces. Consider the case where an interface 
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is defined as a function of x as z(x). The standard way 
of looking at scaling properties is to calculate the inter- 
face width jl7) or standard deviation, i.e. the so-called 
root-mean-square (RMS) roughness, 

where z denotes the mean position of an interface and 
( ) the disorder-average over the different random con- 
figurations. If the interface is self-affine, w should 
scale with L^, ( being the roughness scaling expo- 
nent. For self-affine interfaces the scaling exponent £ 
is expected to be valid also for higher-order statistics. 
This is seen for the height-height correlation functions 
Gk(l) — (\z(x)z(x + l)\ k ) and is applicable to the local 
width as well. Note that there is no a priori reason to 
use a Family- Vicsek-type of scaling Ansatz with a corre- 
lation length |l8) , since there is no dynamical length scale 
for these interfaces. This is an assumption for brittle fail- 
ure, and it will be shown to hold by our data below, and 
is moreover exactly true for perfect plasticity. It is also 
theoretically appealing since slow, adiabatic failure does 
not involve any time scale. 

The local width in two-dimensions is defined analo- 
gously to the global interface roughness w with 

wl c (l) = (j^[z(x) - ZiA, (3) 

where the local interface height zi is averaged over win- 
dows of size I < L. One should note the obvious con- 
nection to G2 that exists both for the local and global 
definitions of interface width. 

The advantage of using more complicated indicators of 
scaling is that one can draw conclusions based on data 
in a much more limited system size range than with the 
global interface width. Of course, such a posteriori tech- 
niques are most commonly used in the context of char- 
acterizing experimental fracture surfaces. Here we note 
the fact that for small L finite size effects make it rather 
difficult to determine the roughness exponent (if one as- 
sumes the interfaces to be truly self-affine to begin with) . 
This is especially true for 3D systems for which the com- 
putational costs become prohibitive easily. 

The global interface roughness w as a function of sys- 
tem size is compared between directed polymers and brit- 
tle fracture interfaces with the dilution-type disorder and 
uniform distribution of J c in Fig. |^. As expected, for 
small L the systems suffer from finite size effects hav- 
ing exponent greater than £ = 2/3 but eventually the 
exponent becomes comparable to the value one obtains 
by fitting a power-law to the large- L-data. Specifically, 
in the dilution case, Fig. ||(a), the effective exponent for 
the fracture surfaces is £fs ~ 0.82 and for the mini- 
mum energy interfaces with exactly the same random 



threshold configurations Cme.< — 0.74. For larger sys- 
tem sizes, which we are able to study numerically only 
in the plasticity limit with large enough number of real- 
izations, Cme,> — 0.67. For the fuses from the uniform 
random distribution, Fig. ^|(b), the fracture surfaces scale 
with (i?s ~ 0.73 and the yield surfaces from the net- 
works with exactly the same random configurations have 
Cme,< — 0.74. Hence the finite sizes effects seem to be 
more similar between the processes than in the dilution 
case. For the larger system sizes of minimum energy in- 
terfaces Cme,> — 0.69. 

Fig. |] compares the local width of directed polymers 
to the brittle fracture interface for the dilution-type dis- 
order with p — 0.8. For the directed polymers one sees 
that the ( = 2/3 scaling is valid for larger system sizes in 
the region, where the window size I < 1/5L. With open 
boundary conditions the scaling region would be larger. 
However, for smaller system size L? = 100 2 there is a vis- 
ible amplitude difference compared to the larger system 
sizes. The fracture surfaces show similar behavior but 
have larger finite size effects when compared to the yield 
surfaces, and they have a larger amplitude, too, than 
minimum energy surfaces in local and in global width 
scaling. 

Our result supports the conclusion of Ref. Q that 2D 
brittle fracture surfaces are in the directed polymer - 
universality class (£ = 2/3), although due to the stronger 
finite size effects the asymptotic region is harder to reach 
than for yield surfaces. 

B. Roughness statistics 

Next we address the higher-order statistics of fracture 
surfaces. For directed polymers one knows that the end- 
point transverse deviation distribution z(x), (= x(t) in 
the ordinary DP notation), is roughly Gaussian for the 
standard case of one fixed, one free end (see e.g. ||) and 
follows a scaling form P[z(L)] ~ f(z/L^). One can like- 
wise write down a scaling form for the interface energy. 
Next we assume that the brittle fracture interfaces obey 
similar self-affine scaling and study the roughness prob- 
ability distribution P(w, L) as a function of L (we con- 
centrate on the p = 0.8 case of the previous subsection). 
Fig. shows a typical example of such distribution: it 
is not centered around zero and is reminiscent of a log- 
normal or Poissonian distribution. This can be under- 
stood qualitatively since the roughness w is also a mea- 
sure of the non-zero maximum transverse displacement z. 
The figure includes a distribution for ME surfaces from 
the same systems, too. 

As one could see already in the previous subsection, 
the 2D fracture surfaces are rougher than the minimum 
energy surfaces (i.e. assuming self-affine scaling in both 
cases, the amplitude of the roughness wjL> is larger for 
FS). This is visible also here since the distribution of 
P(w) is not only wider but extends to higher values for 
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the fracture case. While assuming that P(w) follows the 
same self-affine scaling law as the P[z(L)] we may study 
the disorder-standard deviation 

a(w) = ((w-w) 2 ) 1/2 , (4) 

where w is from a single random system and id is 
the global disorder-averaged roughness calculated using 
Eq. (U). a(w) scales with L 2 / 3 , too, although the data is 
more scattered due to the fact that higher-order statis- 
tics are always more vulnerable to the finiteness of the 
statistics than the integral of them. 

In Fig. |we collapse the data of the cumulative sums of 
the distributions P(w, L) for various L. For both kinds 
of interfaces the data-collapse is better with exponents 
C > 2/3, which is due the finite size effects. For frac- 
ture surfaces £ = 0.8, as in Fig. 0(a), seems to work 
nicely, like for the yield surfaces £ = 0.74 would be bet- 
ter than £ = 2/3. In our opinion the figure justifies the 
assumption of an asymptotically self-similar scaling of 
P[w(L)] ~ f(w/Li). 

C. Scaling of overlap quantities 

The average overlap Po of fracture and minimum en- 
ergy surfaces as a function of the system size for the di- 
lution case is seen in Fig. 0(a). Overlap is defined as 
the fraction of the disorder realizations in which at least 
one (x, z) coordinate pair is common with the fracture 
and yield surfaces. Clearly the overlap is reduced as a 
function of the system size. This is not surprising, be- 
cause with increasing system size the probability of the 
first breaks taking place at the globally weakest place de- 
creases. However, if one expects the overlap quantity to 
originate from the result of depositing the surfaces ran- 
domly (as particles of finite width on a ID line segment of 
length L), one obtains P ran = (A fs wfs+Ameu>me) I 'L- 
Aps and Ame are prefactors needed to compute the typ- 
ical geometrical extent from the roughness values wfs, 
wme- In the figure we have plotted P ran from the same 
RFN configurations as Po, with Aps = Ame = 7.5, 
which is a rather large value to be realistic, hence Po ^ 
P ran - In the figure P ran has a slope ~ —0.2 while the 
asymptotic scaling according to the random deposition 
argument should be P ran ~ L^ 1 ^ 3 . For small systems 
the average overlap is very large, of the order of 0.5 in 
the particular case studied here. 

Fig. 0(b) shows the average size (s) of overlapping clus- 
ters. The overlapping cluster is defined as the number of 
the neighboring common (x, z) coordinate pairs. The 
overlap cluster size saturates at (s) ~ 8.5. Fig. 0(c) 
shows the average total length of the overlap in con- 
figurations, which do have an overlap, i.e. (J2 s )/No, 
where No — PoN . One may write (J2 s )/No = 
((s}N s No)/No = {s}N s , where N s is the number of over- 
lap clusters in a system, which has overlaps. By assum- 
ing that N s ~ L, and (s) =const, we get 8) /No ~ L, 



which is demonstrated in the figure. Since in the large L 
limit s)/No = C1L + C2, C% ~ 23, there is a crossover 
in the systems of size L ~ 27, because in the small system 
size limit ^2 s niust be smaller than L. On the other hand 
C\ = 0.15 tells us that approximately 15% of the length 
of the fracture and yield surfaces are overlapping with 
each other. The scenario is that if the fracture happens 
to start from the same minimum energy valley where the 
DP is located, it will naturally stay localized there; how- 
ever, the associated surface stiffness is weaker and thus 
the excursions. Notice the saturation of the cluster size, 
which agrees with the scenario. 

Fig. shows four examples of what happens with vary- 
ing disorder strength SJ/Jo for L = 100. The subplots 
demonstrate several effects. For the weakest disorder, 
both the interfaces are nevertheless "rough" (i.e. not 
flat), which shows that in spite of a single, growing crack 
even the brittle fracture case can produce a crack which 
fluctuates in the transverse direction. The qualitative 
behavior is the same for both the cases, note that for 
the yield surfaces one expects a Larkin-lengthscale on 
which the interfaces look flat due to the competition be- 
tween disorder and elasticity. With increasing disorder 
the crack is finally localized in the lower part of the sys- 
tem - the threshold field is rescaled in all the cases with 
the "initial" random number being kept constant. Mean- 
while the damage for brittle fracture grows strongly. No- 
tice how the yield (minimum energy) surface moves in 
the system as SJ/J is being changed. For the two cases 
with weakest disorder the surface stays the same. In two 
of the subplots there is considerable overlap between the 
fracture and the yield surfaces: SJ/Jq — 1 leads to a 
total overlap of s — 82. 

Fig. shows the dependencies of the roughness and 
overlap quantities on the disorder strength. Both for the 
dilution case, Fig. 0(a), and for the systems with the 
randomness from the uniform distribution, Fig. 0(b), the 
fracture and yield surfaces are always rough, except for 
the finite size effects in the small 5 J/ Jq limit. Even in 
this case the strong disorder fixed point is attractive and 
we simply have that the system size is smaller than the 
Larkin length above which the asymptotic behavior is 
seen. The roughness increases with decreasing p until 
the bond-percolation limit p c = 1/2 is reached and the 
surfaces become fractals, with the corresponding hull ex- 
ponent. In the insets the average overlap Po and the 
average overlapping cluster size (s) are shown. Po in- 
creases for both type of disorder with increasing disorder 
strength, except in the infinitesimal disorder p = 1 — e 
and 5J/Jq = e limits, where it naturally diverges, the 
same is true for (s) even with finite disorder. In order to 
compare Po with the random deposition argument, P ran 
is plotted from the data of the same configurations with 
Aps = Ame = 8 showing again Po ^ Pran- (s) seems to 
saturate with increasing disorder strength for both types 
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of disorder around (s) ~ 8 — 10, which might be a co- 
incidence, since one could guess it to be disorder-type 
dependent. 



V. DISCUSSION 

This paper has explored the connections between brit- 
tle fracture and minimum energy surfaces. We have given 
numerical support for the idea of these being in the same 
universality class in two spatial dimensions. This argu- 
ment is based on the scaling of interface width, local 
roughness and the statistics of ensembles of interfaces. 
For both scalar brittle fracture and perfect plasticity, or 
minimum energy, interfaces these turn out to have similar 
scaling properties, indicating that the brittle fracture in- 
terfaces have a roughness exponent of 2/3 (the directed 
polymer one) and are also truly self-affine. In spite of 
the fact that we have studied only two types of disorder 
distributions, we nevertheless believe that the numerics 
points to a picture of asymptotically rough cracks in spite 
of the stress-enhancement effects, that one would expect 
to play a role in brittle fracture with a large, dominating 
crack. Such is the case in particular for the dilution- 
type, relatively weak disorder, a main part of our study. 
Notice that for weak disorder even minimum energy sur- 
faces tend to be relatively flat as the amplitude of the 
roughness is small. 

The results presented earlier were obtained for systems 
that were governed by extreme scaling- type arguments. 
The role of the disorder can be tested in another way, 
more relevant to standard fracture surface experiments, 
by introducing a notch, or a row of pre-failed fuses. The 
current distribution around the crack tip is a combination 
of the enhancement effect of the crack plus the additional 
fluctuations created by possible off-path or fracture pro- 
cess zone damage. Born-model simulations by Caldarelli 
et al. [|l9) show that self-affine cracks can avoid surface 
tension effects, i.e. they do not straighten out, if started 
from point seeds. In our case the question becomes if 
the effective surface tension of the notch plus the grown 
crack wins over disorder, remembering that the stress 
enhancement for a symmetric crack is largest on-axis. 
This is a necessary mechanism for any self-affine behav- 
ior, whether of the minimum energy surface universality 
class or not. Fig. |l0| shows the interface roughness of yield 
and brittle fracture surfaces for a fixed system size and 
varying notch length. It transpires that for both ME and 
brittle fracture surfaces the notch effect does not equal to 
flattening. Note the earlier arguments that yield surfaces 
have a higher surface tension. This is again due to the 
memory effect, which renormalizes the thresholds in the 
tangent problems so that they are the smallest on the 
crack axis. 

There are several experimental indications of the con- 
clusion that brittle fracture interfaces exhibit self-affine 



behavior, with a roughness exponent close if not equal 
to the perfect plasticity one. Experiments done on real 
materials can bridge the gap between the two extreme 
limits. For such studies the expected behavior would in 
2D, assuming slow failure, be self-affine as well. The ex- 
traction of the roughness exponent has here been done 
using the local width as a measure. For ensembles of ex- 
periments one should note the statistical implications of 
the scaling of the roughness distribution width and of the 
shape of the width probability distribution. The relative 
"irrelevance" of a notch hints about the possibility of pin- 
ning center -like scaling properties as should be true for 
the perfect plasticity case: the notch pins the final crack 
with certainty if it is large enough. 

Finally we note that there is no rigorous theoretical 
argument that would explain why brittle fracture seems 
to follow ME-type scaling. Indeed we have here studied 
only the scaling of the interface roughness, and the as- 
pect of interface energetics in terms of e.g. the energy 
fluctuation exponent 6 has been left aside. Notice that 
the bare strength properties are governed by logarith- 
mic effects in the case of the brittle failure. For brittle 
fracture interfaces to result from global optimization the 
initial failure thresholds have to be renormalized by the 
correlations which the stress-intensities of the crack his- 
tory induces. For such extremal statistics processes no 
theory exists for the time being, unlike for the case of 
the quenched Laplacian breakdown model for which one 
can use real-space renormalization [ pp[ . 
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system size (L) 

FIG. 2. Scaling of energy of ME surfaces, closed diamonds, 
and fracture quantities: total damage, i.e. number of bro- 
ken fuses nb, open circles; breaking voltage Vb of the network, 
open squares; and breaking current lb of the network, open 
diamonds, as a function of the system size L. The disorder is 
dilution type with p = 0.8. The number of realizations N for 
the brittle failure case is shown in Table |. For ME surfaces 
N = 5000 for L 2 = 10 2 — 50 2 , N = 1000 for L 2 = 60 2 - 350 2 . 
The lines are least squares fits, linear for E and Ub\ L/\/\nL 
for Vb and lb, to the data. 



♦ Fracture surface 
o Minimum energy 




FIG. 1. The voltage-current diagram of a fuse or a me- 
dia. The < > arrow describes an ideal elastic behavior. The 

dashed line describes an ideal elastic break down at the criti- 
cal current ib with the corresponding voltage Vb. The dashed 
arrow describes an perfect plasticity with yield current i y , and 
the dotted arrow describes an elastic-plastic behavior with an 
irreversible yield strain Vi r — Vb at Vi T . The current as a func- 
tion of the voltage with the corresponding elastic or plastic 
behavior of the media may only increase or decrease in the 
directions noted by the arrows. 
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FIG. 3. The interface width w versus the system size L 
for the brittle fracture interfaces, closed diamonds, and the 
minimum energy ones, open diamonds, from the same ran- 
dom networks. Open circles are minimum energy interfaces 
from larger system sizes. The disorder is (a) dilution-type 
with p — 0.8 and (b) from uniform distribution of fuse 
thresholds with S J /Jo = 1. The number of different ran- 
dom realizations for the exactly same fracture and yield sur- 
faces is shown in Table |. ME surfaces have iV=1000 for 
L 2 = 300 2 ,400 2 , JV=500 for L 2 = 500 2 - 700 2 , and iV=200 
for L 2 = 800 2 - 1000 2 in (a), 7V=500 for L 2 = 120 2 - 400 2 , 
and iV=200 for L 2 = 500 2 - 1000 2 in (b). The lines are guides 
to the eye with a slope £ = 2/3. 
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FIG. 4. The local interface roughness wi oc versus the win- 
dow size I for the dilution type of disorder with p — 0.8. The 
system size L 2 = 1000 2 has N = 2000, while for all the other 
system sizes the data is from the same configurations as the 
data in Fig. ||(a). The line is a guide to the eye with a slope 
£ = 2/3. For the minimum energy interfaces the £ = 2/3 
scaling is seen in a region I < 1/5L for larger system sizes, 
while L 2 = 100 2 has a visible amplitude difference. The pe- 
riodic boundaries are used in the transverse direction of the 
external voltage. The correlation between the local width of 
the elastic fracture surfaces and plastic yield ones is clearly 
seen. 
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FIG. 5. The histogram P(w) of the roughness for the ideal 
elastic and perfect plastic yield surfaces in systems of size 
L 2 — 350 2 and dilution type of disorder with p — 0.8. The 
data is from the same configurations as the data in Fig. ^|(a). 
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FIG. 6. Cumulative sums of P(w) for both fracture (a) and 
minimum energy surfaces (b) for various system sizes with di- 
lution type of disorder, p — 0.8. The data is from the same 
configurations as the data in Fig. fcaa). The data has been 
scaled with L 0,8 in (a) and with L 2 ' 6 in (b). 
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FIG. 7. (a) Fraction of the disorder configurations, Po, in 
which the fracture surface and the minimum energy interface 
do have an overlap, i.e. at least one common (a;, z)-coordinate 
pair, open circles. The disorder is dilution type with p — 0.8. 
The data is from the same configurations as the data in 
Fig. ^(a). The closed circles are for the comparison the value 
Pvan = 7.5(wfs + Wmb)/L from the figure ||(a), see the text 
for details, (b) The average size (s) , i.e. the number of com- 
mon neighboring (a;, z)-coordinate pairs, of the overlapping 
clusters as a function of the system size. The overlap cluster 
size saturates at (s) ~ 8.5. (c) The total length of overlap in 
configurations, which do have an overlap in their interfaces, 
si No- The lines are linear least squares fits to data. 
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FIG. 9. Interface width w of the fracture and minimum 
energy surfaces with varying the disorder strength for the di- 
lution type of disorder (a) and uniform distribution of fuse 
thresholds (b). The system size L 2 — 100 2 for each system 
and the number of realizations is shown in Table |. The in- 
sets show of the same systems the fraction of the overlap- 
ping disorder configurations Po, open triangulars, P ra n with 
Afs = Ame — 8 in both cases, closed triangulars, and the 
average overlapping cluster sizes (s) , diamonds, as a function 
of the disorder strength. 
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FIG. 10. Roughness of fracture and minimum energy sur- 
faces in a notched sample with varying the notch length. The 
system sizes L 2 = 100 2 and the disorder is dilution type with 
p = 0.8. 



TABLE I. The number of realizations N performed in sim- 
ulations for exactly the same randomness of brittle failure and 
plasticity. 
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